Short-range correlations in dilute atomic Fermi gases with spin-orbit coupling 
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We study the short-range correlation strength of three dimensional spin half dilute atomic Fermi 
gases with spin-orbit couphng. The interatomic interaction is modeled by the contact pseudopoten- 
tial. In the high temperature limit, we derive the expression for the second order virial expansion of 
the thermodynamic potential via the ladder diagrams. We further evaluate the second order virial 
expansion in the limit that the spin-orbit coupling constants are small, and find that the correlation 
strength between the fermions increases as the forth power of the spin-orbit coupling constants. At 
zero temperature, we consider the cases in which there are symmetric spin-orbit couplings in two or 
three directions. In such cases, there is always a two-body bound state of zero net momentum. In 
the limit that the average interparticle distance is much larger than the dimension of the two-body 
bound state, the system primarily consists of condensed bosonic molecules that fermions pair to 
form; we find that the correlation strength also becomes bigger compared to that in the absence of 
spin-orbit coupling. Our results indicate that generic spin-orbit coupling enhances the short-range 
correlations of the Fermi gases. Measurement of such enhancement by photoassociation experiment 
is also discussed. 

PACS numbers: 05.30.Fk, 05.70.Cc, 67.10.-j, 67.85.Lm 



I. INTRODUCTION 

Recent experimental advances in generating synthetic 
gauge fields are motivated by simulating chargedparti- 
cles in solid state systems by neutral atoms [l|, Q. In 
the presence of external magnetic fields, the degener- 
acy of a manifold of the atom's hypcrfine spin states is 
lifted. The coupling between the manifold of the hy- 
pcrfine spin states and external laser fields gives rise to 
dressed states. The adiabatic elimination of the high en- 
ergy dressed states results in a low energy effective Hamil- 
tonian in which synthetic gauge fields emerge. By such 
schemes, uniform vector potentials [3|, synthetic mag- 
netic Q and electric fields [5| are realized in condensates 
of ®^Rb atoms. With the magnetic field's strength and 
the laser frequency fine tuned, a spin-orbit coupling bilin- 
ear in momentum and pseudo-spin operator components 
in one direction is engineered as well Q. The possibil- 
ity of inducing spin-orbit couplings in two |7| and three 
directions [g is further discussed theoretically. Similar 
attempts to synthesize gauge fields in atomic Fermi gases 
are under active experimental exploration [9|. 

It is an interesting question how the introduction of 
spin-orbit coupling would affect the correlations of di- 
lute atomic gases. In the BEC-BCS crossover problem. 
Tan noticed that the correlations in a homogeneous di- 
lute two-component atomic Fermi gas have an asymptotic 
form [l3| 



(4(r)^l(0)Vi(0)^t(r)> =c(--- 



(1) 



in the regime where r is much less than d the mean dis- 
tance between the particles and bigger than tq the range 
of the interatomic interaction potential U{r). Here f/'o- 
are the field operators for fermions and a^ is the s-wave 
scattering length. The correlation (or contact) strength 



at short distance C has been shown to be linked with 
thermodynamic quantities through a series of remarkable 
relations Il(]|-[l2i. named as Tan's relations, one of which 
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where m is the mass of particles and the free energy den- 
sity is / = — riog(Tre~'^/"^)/y with T the temperature 
and V the volume of the system. 

In the presence of spin-orbit coupling, Eqs. ([l]) and 
([2|) should hold when ro is much smaller than the length 
scale corresponding to the spin-orbit coupling strength. 
For specific, let us consider the spin-orbit coupling of 
the form hso = X^i=x v z '^iPi'^i^ where ai arc the Pauli 
matrices and Kj are the spin-orbit coupling constants. 
Originally, in the absence of hso, Eq- ([I]) can be de- 
rived from the observation 1^ that given the scale sep- 
aration To <C rf in dilute Fermi gases, when one writes 
(V'|(r)?/'](0)i/';(0)V't(r)> = Cx^ir), in the regime r < tq, 
x(r) satisfies the Schrodinger equation of the relative mo- 
tion of two interacting fermions with zero energy 
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The asymptotic form x(^) ^ (!/'' ^ ^/o-s) for r > tq 
is required to connect with the behavior of xi^) in the 
regime r < tq which is solely determined by U{r); as pa- 
rameterizes the effects of U{r) on x(r). Note that since 
we assume that U(r) is not fine tuned close to any reso- 
nance other than in the s-wave channel, the non s-wave 
parts of (V'l(r)7/'!(O)V'^(O)'0^(r)) is neglected for r < rg 
due to the strong suppression by the centrifugal potential 
barrier. The introduction of hso would modify Eq. ([3|). 
Since experimental values of Ki are about 1/d [a, |9|, the 
inverse of the mean interparticle distance, we expect that 



the correction to x due to hgo is of order Kjro, which is 
negligible in dilute gases. A two-body calculation us- 
ing a square well model potential for equal spin-orbit 
couplings in three directions agrees with our expectation 
[13|. Given x unchanged in the lowest order of ftirg, in 
the same way as used in Refs. [13,11J1j O'^s can show that 
Eq. ^ stands. However, how spin-orbit coupling would 
affect the magnitude of the correlation strength at short 
distance is the problem that we are going to study below. 
In this paper, we consider three dimensional spin 
half Fermi gases with spin-orbit coupling hg^ = 
^i=x V z >^iPi^i- The Hamiltonian of the system is 



H —Ho + Hint + Hso, 

Ha= d^rV*^(r)V*(r)/2m, 



H,nt =9 J dh4'\{r)^l{r)M^)M^), 
H,o^ fd^v^Hr)hso'i'{v). (4) 

Here '^ = {ip^, 4'iY' ■ The bare coupling constant g for the 
contact pseudopotential C/(r) = 5<5(r) defined with the 
momentum cutoff A is related to the s-wave scattering 
length Cs via the renormalization 
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We take h = 1 throughout. In the high temperature 
limit, we derive the second order virial expansion of the 
thermodynamic potential of the system by the ladder di- 
agrams. We further evaluate the second order virial ex- 
pansion perturbatively for small Hi and find that the cor- 
relation strength C increases as the forth power of k^. At 
zero temperature, we consider two cases, equal spin-orbit 
couplings in two or three dimensions, in which there al- 
ways exits a two-body bound state with zero net momen- 
tum no matter the value of Cg [l5| . In the limit that the 
mean interparticle distance is much larger than the size 
of the two-body bound state, by the fact that the leading 
contribution to the ground state energy comes from the 
binding energy of the bound state which fermions pair 
into, we show that C becomes bigger compared to that 
in the absence of spin-orbit coupling. Our results indicate 
that nonzero spin-orbit coupling gcnerically enhances the 
correlation strength of the Fermi gases. Such enhance- 
ment can be detected in photoassociation experiment. 



II. VIRIAL EXPANSION AND LADDER 
DIAGRAMS 

When the temperature T is high, the grand canoni- 
cal partition function, Z = Tre^'^"''^)/'^ with /x the 
chemical potential, can be approximated by a virial ex- 
pansion in terms of the fugacity rj = e^'"^ which is a small 
number. The effects of pairwise interactions first appear 



in the second order virial expansion. For two component 
Fermi gases interacting through a short range central po- 
tential without spin-orbital coupling, the second virial 
coefficient due to interactions has been derived UM 
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where En are the binding energies of two-body bound 
states and bg are the phase shifts in the £th partial waves. 
In the presence of spin-orbital coupling, h^o couples scat- 
terings in different partial waves to each other; this cou- 
pling renders classification in terms of angular momen- 
tum as in Eq. ([5]) impossible. However, since the ladder 
diagrams exhaust the two-body scattering processes, we 
use them to calculate the second order virial expansion 

It is instructive to demonstrate how the ladder dia- 
grams reproduce a result agreeable with Eq. ([5]) for a 
Fermi gas whose Hamiltonian is Hq -\- Hint- The varia- 
tion of the thermodynamic potential 51 = —TlogZ due 
to the ladder diagrams is [l^l 



6n^Y. 



with 

n(P,C) = 

Here ^q = 
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C — Cp/2+q — ?P/2-q 



(8) 



Eq — /i, eq = q^ /2m, fs and / are the Bose 
and Fermi distribution functions respectively. The renor- 
malization Eq. ([5]) has been used to obtain Eq. ([T]). The 
branch cut of the logarithmic function lies on the positive 
real axis of its argument. The contour C wraps the real 
axis of the integral variable C,. 

To order of rf , we neglect the Fermi distribution func- 
tions in H, since they contribute at least an extra fac- 
tor of rj; the argument of the logarithmic function in 
Eq. ([7]) becomes the inverse of the T-matrix in vacuum 
t'^iC- P^/im + 2fi), where 
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After changing the variable (' = C,—P^/Am + 2ii, we have 

^" = E/£/^(C' + ^V4m + 2^)log(i"i(C')). 

(10) 

It is generically true that the singularities of 
log(i~^(P, C,')) on the real axis of C,' are left bounded. We 
deform the contour C to wrap the part of the real axis 
of Q' right to the most left singularity. On this contour 



C, we expand fsiC + P"^ /^ra + 2/i) in the integrand to 
the lowest order of 77, which is of order 77^, as 



Direct evaluation of the above equation gives 
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of (5fi consists of only even orders of mniX. To the forth 
order of k^ 

£--n,.„(P,c' + PV4».-2.)^*-(C')-4/|i 



(C - gV"i)^ 



(C - gVm)5 



(16) 



For the second order virial expansion Sfl to second or- 
der of K, , 



b2{X/as) 



Jo TT 1 + (fca^)^ 
-[1 + Erf(A/\/2^a,)]e^'/2^°'. 



(13) 



Here the thermal wavelength is A = ^J2t: jvnT . Equation 
([T3I1 has been obtained by calculating the partition func- 
tion from the two-body eigenenergies in Refs. [ij, [l9| . 
Given that the contact pseudopotential g5(v) scatters 
only the s-wave, and cot5s(fc) = —^/kag, and ~l/m,a1 is 
the binding energy for the only bound state when ag > 0, 
Eq. (fT3|) agrees with Eq. (|6]). For fixed density n ~ N/V 
with N the total number of fermions, Eq. ([T^ is 



(5f]W = -2-^'WTVX'b2{\/ag) 



(14) 



since n = 277/ A^ to order of 77, 

In the presence of Hso, SO, retains the form of Eq. ([7]) 
with n replaced by 



n.o(P,C) 
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with £;" = ^q 



■'qj 



^q = VE7«^' and 0p,q = 



E,«:f(^'/4 - g2)]/Ap/2+qAp/2_q. The index a = ±1 
picks up different helicity branches of noninteracting par- 
ticles [2[- Note that Fermi distribution functions have 
been neglected in II^o for the same reason as stated be- 
fore. 



III. PERTURBATION IN HIGH 
TEMPERATURE LIMIT 

While to obtain the behavior of the second order virial 
expansion 5^ as a function of arbitrary values of ki re- 
quires a full evaluation of the multi-dimensional integral 
(cf. Eqs. dZ]) and ([15])), in the following, we calculate 5^ 
perturbatively in terms of rnXni in the high T limit in 
which A — > 0. Since the integrand of ligo is invariant un- 
der the transformation Ki — >■ — k^, the perturbation series 
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(17) 



Equation (|17p can be reproduced from the diagrams 
shown in Fig. ([T]). The perturbation 6Cl by Hgo comes 
from the ladder diagrams [20] with the vertex Hso at- 
tached. Since the system in the absence of Hgo are in- 
variant under the reflection of momentum p — )■ — p, the 
diagrams with a single vertex ^^ KiPiCJi attached must 
be identically zero. Of second order of k^, similarly, the 
diagrams proportional to KiKj vanish if i 7^ 7; for i = j, 
two KiPiCFi must attach to the same pair of free particle 
propagators as shown in Fig. ([T}. To sec this point, given 
that the three directions are equivalent, let us consider 
attaching KxPxO'x to the ladder diagrams. We choose 
the free particle propagators diagonalizing the z compo- 
nent of the spin operator ct^ . Since the interactions have 
SU{2) symmetry and KxPxX^x flips the spin by unity, dia- 
grams in Fig. ([ij are the only nonzero ones contributing 
to AQ'-^'' within the ladder diagrams. The class of the 
diagrams represented by the left one in Fig. ([IJ gives 
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?P-Hq)(2:-CP-q) 



^P-q)(2-^P-Hq)3 

t(C-PV4"i + 2M), 



(18) 



and the class by the right gives 



sni 



V 



E 



^(Pi/i-q!) 



xt{C-P^/4:m + 2ii). (19) 

Here C. and z are the bosonic and fermionic Matsubara 
frequencies respectively. Of order 7;^, the sum of Eqs. ([THj) 
and HI]) equals Eq. PT]) . 




FIG. 1: The Ladder diagrams contributing to the variation 
of the tliermodynamical potential of the oder of the square of 
the spin-orbit coupling. The point vertices are g, which are 
connected by a pair of free particle propagators of Ho- The 
wiggling vertices correspond to Hso- 



For fixed density n, since n = 2r][l + 
^^(mAKi)^/47r]/A^ to first order of 77 and second 
order of k,, Sn^°^+Sn^^^ = -2-^/^TVX^n%{X/as); the 
second virial coefficient is still &2(A/as), not changed to 
the order of Atf . Consequently, the correlation strength C 
does not change as well. This absence of modification of 
C can be understood in the following way. The contribu- 
tions to 60, proportional to k| from the three directions 
are additive and of the same form. We can reproduce 
the coefficient of k^ in Eq. PT|) by introducing the 
perturbation Hamiltonian Hp = J (fir'^{r)KzPzaz'^{r) 
to Hq + Hint- Note that the unitary tranformation 



explimuz J d^rz{'ipl{r)ip-f{r) 



U ~- 

transforms K — Hq + Hint 

UKU^ =K - {mKl/2)N. We have 
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fiN as 



fl{T, n, Kz) = -Tlog(Tre--^/'^) ^ n{T, n + mKl/2, 0). 

(20) 

For the second order virial expansion 5Vl of order k^ , 



Eq. ([16)) . the non-cross terms ex nj give 
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(22) 



a result expected from expanding Eq. (PH)) in terms of k^ . 
The cross terms oc kJh^ for i y^ j are 
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After changing the integral variable z = 2'!t(' /T and in- 
tegrating by part, we have 




,,2^^rrf^ f d^ _,/2. 
J I A3 A r 27ri^ 
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1 1 2i \ \ ' 



(24) 



6z (a -\' i^fzY 37r ^/z {a + i\fz) 



with a = X/ois- 

Direct evaluation of the contour integral yields 
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where 
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which agrees with Eq. (flT)) . For the correlation strength 
C, let us write its correction of second order Ki as 
A^^(mAKi)^. In the case of with spin-orbit coupling 
only in z direction, 5C = A^ijiXkz)^ ■ However, accord- 
ing to Eq. (|20p , the effect of the spin-orbit coupling only 
in one direction is equivalent to shifting /i. Physically, for 
fixed density, a chemical potential shift should not affect 
C at all; one concludes A = 0. 

The effects of the spin-orbit coupling on C can be re- 
vealed by calculating 6il to the forth power of Ki. From 



Figure ^ shows F{X/as) as a function of X/us, positive 
and well-behaved everywhere. 
To the same order of k, 
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FIG. 2: (Color online) F{X/as) vs A/ss calculted from 
Eq. ^. 
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FIG. 3: (Color online) T{X/as) vs X/gs calculated from 
Eq. (O- 



the thermodynamical potential variation in terms of n is 
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By Eq. ([2) the variation of C is 



(5C = 2-^''^\^n^ 






with 
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(30) 



Figure ([3]) shows that F is always positive; the correla- 
tion strength increases as the forth power of k^ in the 
Hmit mXKi — > 0. Note that since 5C a X]i<7- '*?'*jj the 
correlation strength changes only if spin-orbit coupling 
constants are nonzero at least in two directions. 



IV. ZERO TEMPERATURE 

At zero temperature, we consider that there is a sym- 
metric spin-orbit coupling in cither two or three direc- 
tions, i.e. , Kx = Ky = K and Kz = 0, or Kx = Ky = 
Kz = K. It has been shown that for the two cases there is 
always a two-body bound state with zero center of mass 
momentum for all Os |15| . The eigenenergy of this bound 
state e is given by 



m 
47ras 



n,o(o,e) = o. 



(31) 



where in the expression of IIso the chemical potential is 
set zero. In the limit that the mean distance between 
particles d ^ n~^^^ is much larger than the size of the 



bound state lb, the system can be considered as primar- 
ily consisting of bosonic molecules which two fermions 
pair into. These molecules condense into the zero mo- 
mentum state. Thus we expand in terms of the small 
number £;,n^'^ and have the leading contribution to the 
energy of the system coming from the binding energy of 
the molecules as E = Ne/2. 

For Kx = Ky = K and Kz ~ 0, Eq. ([3T|) reduces to [l^ 
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(32) 



The binding energy Ef, is defined with respect to the 
scattering threshold energy Eth = —rriK^ as Eh = \e 
Eth I . The wavefunction of the bound state is 



V'6(r) 

(^.(r)i n - it) + v-aWi tt) + cwi m , 

(33) 

where ?/'s(r) = - X]k,a cos(k-r)/[2(ek+aAk)-e], '0Q(r) = 

« Ek,a ae^*'^'' sin(k • r)/[2(ek + aAk) - e] with 0k the 
azimuthal angle of k, and the normalization factor M = 
Ek,a2y/[2(ek + «Ak)-e]2. 

It is instructive to check that C extracted from the 
correlation at short distance through Eq. ([T|) satisfies 
Eq. ([2]). By the wavefunction f/'h. 
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= (47r2)-i7V-V,V7V. 
On the other hand, from Eq. pT|) 
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which matches Eq. ([34|) . 

In the case of k^ = Ky ^ k, ^ k, from Eq. pip the 
binding energy E^ is given by [15| 



(36) 




the bound state wavefimction is 
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where the subscript means that f is the spin quan- 
tization axis. The normahzation factor is A/" = 
2n[m'^K'^/{mEbf/^ + l/{mEby/^]. It is also straight- 
forward to show that C extracted from the correlation 
function maintains Eq. ^. 

In the absence of spin-orbit couplings, C = n/A-nas in 
the BEC limit n^/'^Os — >• 0+. The correlation strength 
decreases as as increases and acquires a universal value 
C w 2.7 X (37r2?7,)4/3/367r'' at unitarity l/a, = 0. In the 
BCS limit ri^/^as — >■ 0^, C = alv? /A. For the two cases 
considered above, spin-orbit coupling ensures the exis- 
tence of the two-body bound state of zero net momentum. 
In the limit IbU^'^ — >■ 0, C is contributed primarily from 
the bound state that the fermions pair into, and thus is 
proportional to the density n. In the limit Kas — >■ 0+, for 
symmetric 2D couplings, Eb = l/ma^ + mK^ — m^K'^a^/3 
and C = {n/4nas){l + m'^K'^al/6); for symmetric 3D 
couplings, Eb = 1/ma^ + 2mK^ — m^K^a^ and C = 
{n/ATTas){l+m'^K^as/2). The variation of C proportional 
to K^ and ratio between the coefficients of the k^ terms is 
expected from perturbative calculations as done before at 
high temperatures. Generally, combining Eqs. ([2]), ((32)) 
and ((5S|). we have 



lim 



Tl-i-O C{Kx = Ky ^ Kz ^ k) 
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0, 



< 1, for Gs > 0; 
for flc < 0. 



(38) 



Figure ^ plots 7 = lim„_j.o C{kx = Ky = k^ = 
0)/C{kx = Ky = k^ Kz = 0) versus l/nas- Thus C in- 
creases in the presence of symmetric spin-orbit coupling 
in two or three directions. 



V. DISCUSSION 

Our results at high temperature and at zero tempera- 
ture indicate that generic spin-orbit couplings enhances 
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FIG. 4: (Color online) The ratio 7 vs I/kOs 



the correlation strength C. This enhancement is the net 
result of the two- fold effects brought about by hgo'. the 
change of the density of states of noninteracting particles 
p{io) and the mixing of scattering in different channels. 
The former effect is clearly reflected in determining the 
zero net momentum two-body bound state with symmet- 
ric spin-orbit coupling in two or three directions. For the 
two cases, Eq. ([31]) involves only the scattering between 
two particles of the same helicity and can be written as 
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The change of p(w), especially in the limit w — > from 
~ \/uJ to ^ const, and ~ 1/a/w respectively, gives rise to 
the two-body bound state no matter the value of a^. The 
existence of this bound state guarantees the increase of 
C in the low density limit at zero temperature. In the 
situation where the scattering between two particles with 
nonzero net momentum shall be taken into account, scat- 
terings in intra- and inter-helicity channels are coupled 
together. This mixture can be seen by expressing Hint 
in terms of the fermion operators diagonalizing Ha + Hso 
[2l|. Terms in Hso which determines the second order 
virial expansion correspond to different channels. Re- 
cently Ref. [22] employ the BCS mean field theory to 
study the ground state of attractive spin half fermions 
with spin-orbit coupling, and find that the BCS pairing 
gap A increases with spin-orbit coupling constants. This 
finding agrees with ours since within the BCS theory, 
C = -{m/AT:)[d{E/V)/da-^] - A^. 

The enhancement of C in Fermi gases with spin-orbit 
coupling can be measured by photoassociation experi- 
ment. Previous experiment [23|, across the Feshbach res- 
onance between the lowest two hyperfine states of ^Li 
at magnetic field 834G, associates a pair of Fermi atoms 
in the closed channel of the Feshbach resonance into a 
molecular state. The resultant molecules lose from the 
trap confining the Fermi gas. The loss rate of atoms R 
has been shown to be proportional to C [l2, H, ^3, ■ In 
the presence of spin-orbit coupling, the relation between 



the loss rate R and the correlation strength C should 
remain intact since the photoassociation only involves 
physical processes at a distance ^ tq much shorter than 
the scale introduced by the spin-orbit coupling constants. 
The change of C shall be clearly reflected in the atom loss 
rate. 
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